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Context. The X-ray source 2S01 14+650=LSI+65 010 is a binary system containing a B-type primary and a low mass companion believed to 
be a neutron star. The system has three reported periodicities: the orbital period, P„rb ~11.6 days. X-ray flaring with Pfiare ~2.7 hours and a 
"superorbital" X-ray periodicity P„,pf,- =30.7 days. 
(-H _ Aims. The objective of this paper is to show that the puzzling periodicities in the system may be explained in the context of scenarios in which 
O ^. tidal interactions drive oscillations in the B-supergiant star. 

' Methods. We calculate the solution of the equations of motion for one layer of small surface elements distributed along the equator of the 
star, as they respond to the forces due to gas pressure, centrifugal, coriolis, viscous forces, and the gravitational forces of both stars, which 
provides variability timescales that can be compared with those observed for 2S0 114+650. In addition, we use observational data obtained at 
the Observatorio Astronomico Nacional en San Pedro Martir (OAN/SPM) between 1993-2004 to determine which periodicities may be present 
in the optical region. 

^ ■ Results. The models for circular orbits predict "superorbital" periods while the eccentric orbit models predict strong variations on orbital 
k>( timescales, associated with periastron passage. Both also predict oscillations on timescales of ~2 hrs. We suggest that the tidal oscillations lead 

j_j to a structured stellar wind which, when fed to the neutron star, produces the X-ray modulations. The connection between the stellar oscillations 
and the modulation of the mass ejection may lie in the shear energy dissipation generated by the tangential motions that are produced by the 
tidal effects, particularly in the tidal bulge region. From an observational standpoint, we find indications for variability in the He I 5875 A line 
on ~2 hrs timescale and, possibly, the "superorbital" timescale. However, the line profile variability exceeds that which is predicted by the tidal 
interaction model and can be understood in terms of variable emission that is superposed on the photospheric absorption. This emission appears 
to be associated with the B-supergiant's stellar wind rather than the vicinity of the companion. 

Conclusions. The model calculations lead us to conclude that the B-supergiant may be the origin of the periodicities observed in the X-ray 
data, through a combination of a localized structured wind that is fed to the collapsed object and, possibly, by production of X-ray emission on 
its own surface. This scenario weakens the case for 2S01 14+650 containing a magnetar descendent. 

Key words, interacting binaries - oscillations — stellar rotation — LSI+65 010 
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1. Introduction 

The X-ray source 2S0114+650 (Dower & Kelley 1977) is a 
binary system containing a luminous B-type star and a low 
mass companion believed to be a neutron star. Unlike many 
such systems, however, the short-period X-ray pulsations asso- 
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ciated with a rapidly rotating neutron star are apparently absent. 
Koenigsberger et al. (1983) reported ~15 minute pulsations 
which, however, were not found in later observations except 
those of Ginga (Yamauchi et al. 1990). On the other hand, X- 
ray flaring activity was shown to be periodic (Finley et al. 1992; 
Hall et al. 2000) with Py/a,.^, ~2.78 hrs, a period that in more re- 
cent INTEGRAL data has decreased to ~2.67 hrs (Bonning & 
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Falanga 2005). Finley et al. (1992) suggested that this period 
could actually be due to the pulsar's rotation which, however, 
would make it an extremely slow pulsar. Li & van den Heuvel 
(1999) showed that the neutron star in the 2S01 14+650 system 
could have reduced its rotation rate within the lifetime of the B- 
supergiant companion only if it had a very large (> lO'"* Gauss) 
magnetic field at the time of its birth. Thus, the new-born pulsar 
would have been a magnetar. 

An alternative scenario for explaining the 2.7 hour period- 
icity that was also suggested by Finley et al. involves pulsations 
of the optical component by which a structured stellar wind 
with alternating high and low density regions is produced. The 
density variations of material being accreted by the collapsed 
companion would lead to the observed periodic flaring activity. 
However, it is not evident: a) whether a B-supergiant may os- 
cillate with this short period; nor b) how the stellar oscillations 
can produce the variable mass-loss; nor c) how the ejected mass 
can be distributed in the precisely structured stellar wind that 
is needed for the periodic modulation of the accretion rate onto 
the neutron star. If a plausible scenario involving pulsations of 
the B-star could be constructed, the magneiar hypothesis would 
be weakened. 

A second X-ray variability timescale in 2S01 14+650 is the 
orbital period. Corbet et al. (1999) suggested that this variabil- 
ity may be due to eclipses of the X-ray source, but Mukherjee 
& Paul (2006) find that the change in the ASCA X-ray luminos- 
ity can be ascribed to a local change in the density of the stellar 
wind from which the neutron star accretes. The latter would 
occur either if the orbit is eccentric or if the mass-loss charac- 
teristics of the B-star are orbital-phase dependent. Interestingly, 
the description given by den Hartog et al. (2006) of the broad 
maximum in the RXTE-ASM and INTEGRAL data as consist- 
ing mostly of flares supports the idea of a non-steady mass-loss 
rate from the B-star, with larger amplitude variability at phases 
when the average X-ray flux is maximum. This behavior is one 
that might be expected in eccentric orbits. On the other hand, 
however, a third periodic X-ray variability timescale has been 
recently reported by Farrell et al. (2006): Psuper ~30.7 days, 
which they suggest may be caused by a warped accretion disk 
around the collapsed object. However, the absence of He II 
4686 A emission (Reig et al. 1996; Koenigsberger et al. 2003) 
argues against the presence of a significant accretion disk. 

Non-synchronously rotating stars in binary systems are 
well-known for exhibiting numerous periodicities due to the os- 
cillations that are driven by the tidal forces (Press & Teukolsky 
1977; Zahn 1977; Savonije et al. 1995; Kumar et al. 1995). In 
particular, Moreno et al. (2005) found "superorbital" periodici- 
ties in models of binary systems with circular orbits, in addition 
to shorter period variability. Thus, the question naturally arises 
of whether we can produce a consistent picture for explaining 
the puzzUng periodicities of 2S01 14+650 through scenarios in- 
volving tidal interactions. 

In Section 2 we present the observational data; in Section 
3 we summarize the tidal interaction model that we use and 
describe its predictions for an eccentric and a circular orbit; in 
Section 4 we describe the behavior of the He 1 5875 A fine; and 
Section 5 contains our conclusions. 



2. Observations and data reduction 

Nine different epochs of observations of LSI+65 010 are an- 
alyzed in this paper. The data of four of these were described 
in Koenigsberger et al. (2003) three of which (1993 Oct 9, 10, 
13; 2001 Jan; and 2001 Oct) were obtained wifli flie REOSC 
echelle spectrograph on the 2.1 m telescope of the Observatorio 
Astronomico Nacional en San Pedro Martir (OAN/SPM) and 
one was obtained in 1995 with the coude spectrograph on the 
University of Hawaii 2.2 m telescope. The new sets of data 
were obtained with the 2.1m OAN/SPM telescope as follows: 
on UT 2003 Jan 18, 2004 Oct 29-3 1, Nov 19-26, and Dec 15-19 
with the MES-SPM spectrograph, and on 2004 Nov. 2 and Nov. 
5-7 with the REOSC echelle spectrograph. The total number of 
spectra that were analyzed is 182. 

The REOSC echelle camera system is described by Diego 
& Echevarria (1994). With the 300 grooves mm"' echellete, the 
system gives a resolution of of 7.8 A mm ' (0. 15 A pixeL') at 
H j8 and 10.7 A mm"' (0.20 A pixel"') at Ua. The resolution 
is R=17,000, which corresponds to 17 km s"' per 2 pixels. The 
sht was set at a width of 150 jum, which corresponds to 2". 
The camera of the MES-SPM spectrograph is of a slow focal 
ratio, f/7.5, producing a stability in the wavelength dispersion 
of +0.03 A or a velocity dispersion of ±1.5 km s ' (Richer, 
2003). LSI+65 010 was observed tiirough tiie He I 5875 A fil- 
ter, which covers a bandwidth of approximately /1/15860-5903 
A allowing the use of the interstellar Na I D lines at AA5SS9, 
5895 A as absolute velocity standards. Spectra of LSI+65 010 
were interleaved with spectra of the Th-Ar lamp. 

Data reduction was performed with IRAF' and MIDAS^ 
using standard reduction procedures. Table 1 summarizes the 
general characteristics of the data sets. Column 1 identifies the 
epoch of observation, column 2 the instrument, column 3 lists 
the average JuUan date, column 4 the number of spectra, and 
column 5 lists the average signal-to-noise ratio and standard 
deviation about the mean (in parenthesis). This is the average 
S/N ratio of the individual spectra. The standard deviation is 
given only for nights when at least five individual exposures 
were made. 

The MES-SPM instrumental response has a very compli- 
cated shape, due to the interference filter used, with a bump 
very near the red edge of the He I line. If not taken into ac- 
count properly, this bump could be interpreted as a red emis- 
sion component. In order to avoid this problem, the spectrum 
of the radial velocity standard star HD22484 was used to rectify 
the continuum. This was done by cutting out the photospheric 
lines, which are all very narrow in this star, and the resulting 
"cleaned" continuum distribution was used to normahze the 
spectra of LSI+65 010. 

The stability of the wavelength calibration was tested by 
measuring on each individual spectrum at least one of the inter- 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Asociation of Universities 
for Research in Astronomy, Inc.under cooperative agreement with the 
National Science Fundation. 

^ MIDAS is the acronym for the Munich Image Data Analysis 
System which is developed and maintained by the European Southern 
Observatory. 



Koenigsberger et al.: LSI+65 010=280114+650 
Table 1. Radial velocity and residual emission He I 5875 A 



3 



Epoch 


Instnim. 


<MJD> 


N 


S/N 


^orb 


sup 


s.d.(ISM) 

km/s 


RV (±s.d.) 

km/s 


-EW 

A 


19930ct 


REOSC 


49269.54 


5 


26(8) 


0.804 


0.78 


1.1 


-41.5 (1.8) 


0.35 






49270.33 


5 


29(7) 


0.881 


0.80 


1.5 


-37.7 (2.3) 


0.39 






49273.39 


14 


20(7) 


0.145 


0.90 


2.6 


-41.7 (3.9) 


0.45 


1995NOV 


UH coude 


50049.38 


13 


36(5) 


0.052 


0.11 


0.7 


-65.8 (4.6) 


0.61 






50050.34 


9 


33(9) 


0.135 


0.14 


0.4 


-63.6 (3.5) 


0.65 


2001Jan 


REOSC 


51931.18 


11 


84(17) 


0.304 


0.22 


2.1 


-53.9 (3.8) 


0.00 


2001Oct 


REOSC 


52195.36 


12 


68(11) 


0.082 


0.80 


3.4 


-36.9 (4.2) 


1.10 


2003Jan 


MES 


52657.40 


12 


38(5) 


0.920 


0.81 


1.3 


-50.3 (3.3) 


1.09 


2004Oct 


MES 


53307.33 


2 


60 


0.957 


0.92 




-23.3 


0.60 






53308.30 


2 


50 


0.041 


0.95 




-36.5 


0.99 






53309.27 


1 


53 


0.125 


0.98 




-44.0 


0.97 


2004NOV 


REOSC 


53311.34 


15 


60(11) 


0.304 


0.05 


2.4 


-62.3 (5.2) 


0.70 






53314.30 


3 


55 


0.559 


0.14 




-87.6 


0.87 






53315.20 


2 


62 


0.637 


0.17 




-77.4 


0.59 






53316.30 


2 


47 


0.731 


0.21 




-75.0 


0.81 


2004NOV 


MES 


53328.16 


1 


69 


0.754 


0.59 




-47.5 


0.75 






53329.43 


1 


74 


0.864 


0.64 




-42.3 


0.66 






53330.29 


1 


53 


0.938 


0.66 




-38.8 


0.57 






53333.40 


2 


70 


0.206 


0.77 




-63.2 


0.62 






53334.41 


1 


60 


0.293 


0.80 




-70.7 


0.21 






53335.40 


1 


66 


0.378 


0.83 




-80.2 


0.49 


2004Dec 


MES 


53354.23 


16 


58(9) 


0.002 


0.44 


1.5 


-30.1 (1.3) 


0.84 






53355.20 


17 


53(6) 


0.085 


0.47 


1.0 


-43.6 (1.5) 


0.90 






53356.20 


15 


53(13) 


0.172 


0.51 


1.1 


-46.0(1.5) 


0.77 






53357.20 


18 


54(8) 


0.258 


0.54 


1.1 


-59.1 (2.1) 


0.95 






53358.20 


8 


59(9) 


0.344 


0.57 


0.3 


-60.8 (1.6) 


1.15 



stellar lines in the vicinity of He I 5875.80 A. The Na I A5889, 
5895 A ISM lines in LSI+65 010 have several velocity com- 
ponents (Koenigsberger et al. 2003), two of which are narrow 
and clearly resolved in the UH-coude and MES-SPM data. We 
used the highest velocity component of /15889.95 A located at 
~5888.90 A to establish the precision of the velocity measure- 
ments for these spectra. This ISM component is not resolved 
in the echeUe data. Hence, we used the zero-velocity compo- 
nent of Na 1 5889.95 A and the diffuse interstellar band (DIB; 
Herbig 1995) located at /15849.887 A to determine the preci- 
sion of velocity measurements on the REOSC spectra. Column 
8 of Table 1 lists the standard deviation about the average value 
of the chosen ISM features, for all cases in which 5 or more 
spectra were obtained. In general, the precision of the velocity 
measurements in the UH-coude and MES-SPM spectra is bet- 
ter than 1.5 km s \ while the REOSC echelle data allows RV 
measurements with a precision better than ~A km s ^ 

In order to account for possible zero-point shifts in the 
wavelength scale, and at the same time, correct for the motion 
of the Earth and Sun, all radial velocities of He I /15875 A were 
measured with respect to the centroid of the nearby interstel- 
lar absorption lines. In practice, this was accompUshed by: a) 
fixing v=0 km s ' on the velocity scale at either /15849.89 A 
/15888.90 A or ^5889.95 A; b) measuring the velocity of the 
ISM component and the He I 5875 photospheric line; c) sub- 
tracting the measured velocity of the ISM line from the He I 
line velocity. This procedure was followed for each individual 
spectrum. The ISM lines were measured by fitting one or sev- 



eral Gaussian functions while the centroid of the He I 5875 A 
Une was measured with IRAF by integrating the pixel values 
between the two edges of the absorption Une profile that inter- 
sect the level 0.9 of the normalized continuum. This method 
for He I 5875 A was chosen primarily because it reduces the 
effects produced by the variabiUty of the extended wings of the 
line. Column 9 of Table 1 Usts the average value of the RV, 
obtained by averaging the values of the individual RV mea- 
surements of He I corrected for the shift in the ISM line, for 
each night (+s.d. in parentheses). This velocity is with respect 
to the laboratory wavelength, adopted as /lo =5875.80, which 
is the average wavelength of the two tabulated He I lines at 
/1/15875.62 and 5875.97 A. 



3. The tidal interaction model 

Moreno & Koenigsberger (1999) and Moreno et al. (2005) pre- 
sented a simple model that describes the behavior of a star's 
surface layer in the presence of the perturbing effect of a binary 
companion. The model calculates the solution of the equations 
of motion for one layer of small surface elements distributed 
along the equator of the star as they respond to gas pressure, 
centrifugal, coriohs and viscous forces, as well as the gravita- 
tional forces of both stars. The equations of motion are solved 
in the non-inertial reference frame that is centered on Mi, the 
star being modeled, and that rotates with the same angular ve- 
locity as the orbital motion of M2, the companion star. The 
main body of Mi, interior to the thin surface layer, is assumed 
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to behave as a rigid body, with the tidal deformation appearing 
only in the surface layer 

For the numerical calculation, the required input parame- 
ters are Mi and M2, Mi's radius and rotation velocity, Ri and 
Vrof, respectively, the orbital period Pgrb^ the eccentricity e, the 
inclination of the orbital plane /, the argument of periastron 
ojper, the viscosity v, and the relative thickness 5R/Ri of Mi's 
external oscillating layer The output consists of temporal series 
describing, among other characteristics, the maximum pertur- 
bation in stellar radius, R„ma , and radial velocity, Y,„ax, the az- 
imuthal component of the velocity perturbation for all azimuth 
angles at given points in time, AVo', and the photospheric ab- 
sorption line profiles as detected by an observer in an external 
inertial reference system. 

The orbital eccentricity and the stellar rotational velocity 
are two critical parameters for understanding the variability 
patterns in binary systems. If e>0, the perturbation produced 
during periastron passage is a dominant effect and hence, the 
strongest variability occurs on orbital timescales. On the other 
hand, if e=0, "superorbital" or "suborbital" periodicities ap- 
pear. Whether the periodicity is "superorbital" or not depends 
on the ratio of the rotational angular velocity to the orbital an- 
gular velocity, j6o = (jJolQ.per- Here wq -'^J^l^mt =Vro,/Ri is the 
rotation angular velocity of the rigidly rotating inner stellar re- 
gion, and Ri the unperturbed surface stellar radius. Q.per is the 
orbital angular velocity, and its value at periastron is adopted 
for the general case of an eccentric orbit. Thus, y6o is a measure 
of the non-synchronicity between the stellar rotation and the 
orbital angular velocities at periastron, and it can be written as: 



Ri (l-He)i/2 
where Vr„, is in units of km s"' 



(1) 



Ri in units of R0, Po,i in days. 



and e is the eccentricity of the orbit. 
3.1. Application to LSI+65 010 

The spectral line profiles of LSI-i-65 010 have peculiar shapes 
and are variable (Koenigsberger et al. 2003) which results in 
a relatively large uncertainty for v sin /=96+20 km s"' (Reig 
et al. 1996). We illustrate in Figure 1 the theoretical rotation- 
ally broadened line profiles for v sin i -96 km s ' and 120 
km s The cores of both of the observed profiles are indeed 
well matched with the v sin i -96 km s ' theoretical predic- 
tion, but their blue wings are significantly more extended. For 
comparison, the observed line profile obtained by Canalizo et 
al. (1995) of the B-supergiant optical counterpart of Cyg X- 
1, HDE 226868, is also included in this figure. HDE 226868's 
absorption is better matched by the v sin / =120 km s"' theoret- 
ical profile, although the extent of its blue wing is also under- 
predicted. The blue wing of strong photospheric absorptions 
in early-type stars is generally more extended due to the ex- 
pansion of the external atmospheric layers as they accelerate 
and become the base of the stellar winds. Variability in the 
blue wing as displayed in Figure 1 for LSI-i-65 010 indicates 
that the atmospheric properties at the base of the wind undergo 
significant time-dependent changes. These changes occur on 
timescales of ~days (Koenigsberger et al. 2003). 
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Fig. 1. Observed line profiles of LSI-H65 010 at two different 
orbital phases compared with theoretical lines profiles rotation- 
ally broadened to v sin(i)=96 km s"' (dashes) and 120 km s~' 
(dots), and with HDE 226868, the B-supergiant counterpart of 
Cyg X-1. Profiles are shifted by +0.03 on the vertical scale for 
clarity in the figure. 



It is important to note that LSI+65 OlO's profile at =0.304 
shown in Figure 1 substantially resembles HDE 226868's line 
profile, as well as the theoretical profiles. Thus, the =0.304 
spectrum corresponds to a time when the B-supergiant is min- 
imally perturbed and is probably the best profile for estimating 
V sin i . Adopting the value found by Reig et al. (1996) and the 
reported estimates of / ~61-90°, ^rot ~96-109 km s"'. The 
other stellar parameters derived for LSI+65 010 are Mi =16+5 
M0,Ri =37+15 Rq, and M2 =1.7 Mq for the neutron star (Reig 
et al. 1996). 

Keeping Mi =16 Mq, M2 =1.7 Mq and e=0 fixed, we per- 
formed a series of model calculations searching for variabil- 



ity on the P super ~30.7 days, finding that P super increases lin- 
early with ySo for ySo in the range 0.57-0.63. "P super ~30.7 days 
was achieved with jio =0.615. Fixing this value in Eq. (1) con- 
strains Ri to the range 36.2-41.1 Rq, for values of ^roi between 
96-109 km s '. Figure 2 illustrates the characteristics of the 
"superorbital cycle" for two parameter sets (Ri ,Vrn/)=(40 R©, 
106 km s"') and (36.2 Rq, 96 km s"'). A slow increase in the 
size of the tidal bulge over most of the cycle is followed by a 
more abrupt growth to an absolute maximum lasting ~5 days. 
Oscillations on timescales of ~hours also appear and, for the 
cases illustrated in Figure 2, they lie in the range 2-3 hours. 

In contrast to the circular orbit case, there are no "superor- 
bital" periods on <100 day-timescales in eccentric orbit case, 
so the value of ySo is unconstrained except for the limits on v,o, 
and Ri given by Reig et al., which lead to the range ySo ~0.40- 
0.90. An example of the tidal interactions predicted by the 
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Fig. 2. Radial extent of the tidal bulge, R,nax (top), and maxi- 
mum surface velocity, V„mv, plotted as a function of time for 
circular orbit models with /3o =0.615. Results for two different 
stellar radii, Ri =40 Rq (dotted hne; model C-1) and Ri =36.2 
Rq (model C-2) are shown, illustrating the ~31-day "superor- 
bital" period that appears. Superposed oscillations on ~hours 
timescales are indicated by the arrows, though they are barely 
visible on this figure. 



model for ySo =0.66 is shown in Figure 3. The data for one 
orbital cycle are plotted and the abscissa is given in units of 
orbital phase with periastron (0 =0.0) and apastron indicated 
by the vertical tick marks. It is interesting to note that although 
the tidal bulge grows gradually as periastron is approached, it 
reaches maximum extent only after periastron passage due to 
the retarding action of the viscosity. 

The short-timescale oscillations are most clearly seen in the 
bottom plot of Figure 3. Their frequency gradually decreases 
with time and corresponds to periodicities in the range ~ 1.6-3.4 
hours. The predicted amplitude is largest in the phase interval 
~0.15-0.40. 

The idea of tidal oscillations in the tangential direction 
is somewhat less familiar than that of the radial oscillations. 
However, the "swashing" back and forth of the surface ele- 
ments as they approach and pass the sub-stellar point is respon- 
sible for some of the strongest variability that is produced in 
photospheric line profiles. The tangential velocity perturbation, 
A Vo' for the e=0. 16 case is illustrated in Figure 4 for different 
orbital phases, showing that the largest amplitudes occur near 
the primary tidal bulge, O' ~360°. Note also that the perturba- 
tions are strongest after peristron passage, around orbital phase 
0.20. 

In Table 2 we list the observationally-derived parameters of 
LSI-i-65 010 as well as the parameters of the model calculations 



for the eccentric orbit and the two calculations for the circular 
orbit. 

3.2. Viscosity and sliear energy dissipation 

The Ri =40 R© models are highly unstable, which is easy to un- 
derstand considering that this radius is very close to the Roche 
radius, requiring relatively large values of the viscosity to al- 
low the computation to be performed. The oscillation ampli- 
tudes scale inversely with the viscosity, so larger amplitudes 
would be attained if we could use smaller v values. However, 
as described by Moreno et al. (2005), when the oscillation am- 
plitudes become too large, the individual surface elements for 
which the equations of motion are being solved become de- 
tached or begin to overlap, at which time the numerical com- 
putation is halted. To prevent this from happening, the value 
of V may be increased. Our practice is to choose the smallest 
value of V that allows the code to run over at least ~5 orbital 
cycles in order to compute the line profiles after the initial tran- 
sitory phase of the calculation. In the real star, the phenomena 
that we describe as "overlapping or detached surface elements" 
have no significance except that they suggest the possibility of 
large-scale relative motions that may be interpreted in terms of 
the more familiar term macro-turbulence. With this interpreta- 
tion in mind, the value we use for v in the calculations may be 
thought of as describing the conditions of the stellar surface at 
the limit where large-amplitude turbulent motions develop on 
small length scales. 

The viscosity is related to the a-parameter (Shakura & 
Sunyaev 1973) used in studies of accretion disks in Keplerian 
orbits through v ~ (a/cL>)(kT/fj.mH), where T is the tempera- 
ture, yU is the mean molecular weight, and w is a characteristic 
angular velocity. Associating <y with the relative motion of the 
external layer with respect to the rigidly-rotating inner region 
of our model, and assuming T=20,000 K, the values of v we use 
in this paper all correspond to a <0.01, well within the range 
expected for the outer layers of a supergiant star. 

A direct consequence of the behavior of AVo' and of includ- 
ing the viscosity in the calculation is that the surface layer slid- 
ing over the inner, rigidly-rotating body leads to shear energy 
dissipation, E. Note that £ is a potential source of additional 
energy that may be fed into the surface layers of the star and it 
has been suggested that E can contribute towards driving stellar 
winds or mass-ejections, or may play a role in the generation 
of magnetic fields (Koenigsberger et al. 2002; Moreno et al. 
2005; Toledano et al. 2006). All these processes are generally 
associated with X-ray emission. Within this context. Figure 4 
implies a concentration of surface activity on the hemisphere of 
Mi that faces its companion. Hence, in addition to the effects 
produced by the accretion onto the neutron star and the ensu- 
ing X-ray emission, a strong degree of variability due to the 
activity on the B-star surface facing its companion is expected, 
particularly at periastron passage (if the orbit is eccentric). 

The two approximations used in our model calculations that 
must be noted are the artificial division of the star into an outer 
oscillating layer and an inner rigidly rotating body, and the fact 
that the calculation of the oscillations is performed only for the 
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Fig. 3. Radial extent of the tidal bulge, R,„;, i , and maximum sur- 
face velocity, Vma.\, over one orbital cycle for the e=0. 16 eccen- 
tric orbit. The abscissa is orbital phase. Periastron (0 =0.0) and 
apastron are indicated. The short-timescale oscillations, best 
seen in the bottom panel, have periods ~hours that gradually 
increase after periastron. 
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Fig. 5. Radial velocity curves for three different periods: 
Crampton et al.'s (1985) 1 1.591 (top, shifted by -)-120 km s"'), 
den Hartog et al.'s 11.5995 (middle, shifted by +60 km s"') 
and 11.5980 days (bottom). Data plotted are: Crampton et al. 
(1985) (open triangles), Reig et al. (open squares), our 1993 
data (stars) and our 1995-2004 averages (filled squares) of at 
least 5 spectra per night. 
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Fig. 4. Azimuthal velocity perturbation of the surface layer as 
a function of angle for the eccentric orbit case for four different 
orbital phases, as indicated. Note that the largest azimuthal ve- 
locities occur near ^ =0.12, rather than at periastron passage. 
Large values of A V(D' lead to large energy dissipation rates. 



equatorial region of the star The former limits the ability of 
the code to reproduce many of the oscillation frequencies that 
depend on the structure of the entire star, while the latter yields 
upper limits to the oscillation amplitudes since the strongest 
perturbations occur in the orbital plane. Finally, we note that 
the line-profile calculation does not currently include effects 
due to possible variations in the effective temperature over the 
stellar surface. 

4. LSI+65 010: Results from the observational data 

AA. RV curve 

The B-type optical counterpart of 2S01 14+650, LSI+65 
010=V662 Cas, was first identified by Margon & Bradt (1977) 
based on its variable Ha emission and classified by Reig et al. 
(1996) and Liu & Hang (1999) as Bl la. Using measurements 
of photospheric absorption-line radial velocity (RV) variations 
from data obtained between 1978 and 1984, Crampton et al. 
(1985) concluded that the orbit of the B-supergiant was either 
circular, with P"^*^'' =11.591+0.003 days or slightly eccentric 
(e=0.16), with pj'l =1 1.588+0.003 days. Recently, den Haitog 
et al. (2006) found P„,h =1 1.5995+0.0053 days from an analy- 
sis of RXTE-ASM and INTEGRAL data. 

Figure 5 presents a plot (top) of the Crampton et al. (1985) 
data and our average data from Table 1 all folded in phase with 
the circular-orbit period and initial epoch. Our data in this fig- 
ure are limited to points for which 5 or more individual spectra 
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from of a Ml =16 M© primary star tidally perturbed by a 
M2 =1.7 Mq companion in a e =0.16 (dot-dashes) or a cir- 
cular (dashes) orbit. Phase 0.0 is maximum receding velocity 
for the circular orbit and periastron passage for the eccentric 
one. Time of X-ray maximum according to dem Hartog et al. 
ephemeris but with P= 11. 5980 days is indicated. 



were obtained during the night. A significant phase shift be- 
tween the modulation of the Crampton et al. (1985) data and 
the more recent data is evident. This phase shift is reduced if 
the data are folded with the den Hartog et al. P= 1 1.5995 days 
(middle plot in Figure 5). An even better agreement between 
the two data sets is achieved with P= 1 1 .5980 days, which is 
well within the uncertainties given by den Hartog et al. The bot- 
tom plot of Figure 5 also includes the data published by Reig 
et al. (1996), obtained from RV measurements of He I 6678 A. 
They follow in part the general trend of the other two data sets, 
but display significantly larger negative RV at phases ~0. 1-0.4 
which correspond to orbital phases in which the secondary star 
is "in front" of the B-supergiant. This means that the photo- 
spheric absorption lines must have had a significantly more 
extended blue wing during these orbital phases, suggesting a 
difference between the hemisphere of the B-star that faces its 
companion and the opposite hemisphere. Considering the dis- 
cussion of the previous section, we speculate that the difference 
between the two hemispheres may be associated with the tidal 
bulge and its oscillations. 

In Figure 6 we now plot all the individual data measure- 
ments folded with P= 11.5980 days and To =43669.9 MJD. 
This value of To is derived by taking the initial epoch given 
by Crampton et al. (1985) for their circular orbit solution and 



subtracting 40 cycles of 1 1 .5980 days each. To corresponds to 
the time of fastest receding radial velocity, which means that 
the hemisphere of the B-supergiant that faces its companion 
starts coming into view at =0 and is occulted in our line- 
of-sight at (p >0.5. The initial epoch T(,„ of the Crampton et 
al. (1985) eccentric orbit solution defines phase as periastron 
passage. The different initial epochs are all listed in Table 2. We 
indicate with an arrow in Figure 6 the center of the broad X- 
ray maximum reported by den Hartog et al. According to their 
ephemeris, the X-ray "high" state is centered at (p ~0.12 for 
the circular orbit solution or ~0.07 if the orbit is eccentric. It 
is difficult to understand why X-ray maximum should be cen- 
tered at these phases instead of around =0.25 unless the orbit 
is eccentric, in which case X-ray maximum would occur very 
close to the phase for which our tidal model predicts maximum 
perturbation of the tidal bulge (see Figure 2). 

The curves drawn in Figure 6 correspond to the two orbital 
solutions (e=0, e=0. 16) but, compared to the large scatter of the 
data, are nearly indistinguishable. As we will show below, the 
dispersion of the RV data is intrinsic to the line profiles. Thus, 
the RV's contain information additional to the orbital motion 
making the determination of the orbital solution very uncertain. 
It is interesting to note that the RV curves of other HMXRBs 
such as HD 153919/4U1700-37 (Hammerschlag-Hensebergeet 
al. 2003) and HD77581/Vela X-1 (Barziv et al. 2001) also suf- 
fer from large intrinsic scatter, introducing significant uncer- 
tainties in the determination of the masses of their respective 
collapsed companions. 

The individual RV measurements made on the MES-SPM 
data are plotted in Figure 7, showing that there is a systematic 
trend for cyclical variability with a characteristic timescale of 
~2.1 hours. The error bars correspond to the standard devia- 
tion of ISM line measurements, and are significantly smaller 
than the amplitude of the cyclic variation. The similarity be- 
tween this variability timescale and Vfiare is interesting, but we 
unfortunately do not have sufficient orbital phase coverage to 
establish when these cyclical variations are present and when 
they are not. There are two times in Figure 7 when they are ab- 
sent, one of which (0 =0.258), intriguingly, coincides with the 
conjunction corresponding to Mi in back of its low-mass com- 
panion. If the optical ~2-hour oscillations were intrinsic to the 
neutron star region, this would be the most favorable phase to 
observe them. The absence of oscillations at this phase is thus 
perplexing. 

4.2. Line profile variability 

The large scatter in the RV curve of LSI-i-65 010 can be traced 
to line profile variability, as illustrated in Figure 8 where the He 
I 5875 A line profiles are stacked in order of increasing phase. 
The shape of the observed profiles does not repeat with orbital 
phase, nor is there a clear pattern of variability over the or- 
bital cycle. Some spectra display extraordinarily extended blue 
wings (e.g., the 1995 data) and some spectra have an extended 
"red" emission wing (e.g., orbital phases 0.002, 0.344, among 
others). In contrast, the line-profile variability that is predicted 
from the tidal oscillations is periodic with orbital phase, the 
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Fig. 7. Centroids of individual spectra obtained in 2003 and 
2004 with the MES-SPM spectrograph plotted as a function 
of time, with arbitrary initial times. Different nights are given 
a vertical shift for clarity. Error bars correspond to the uncer- 
tainties determined from the ISM-line measurements. The sine 
curves have frequency of 11. 4 cycles day"'. 

strongest variations occuring in the phase interval 0.0-0.3. This 
is illustrated in Figure 9, where we plot the line profiles com- 
puted in the e=0. 16 modeP. The perturbation on the line profile 
first appears at periastron passage on the blue wing and moves 
towards the red wing, becoming imperceptible by orbital phase 
~0.3. This pattern repeats over each orbital cycle. Hence, the 
observed behavior of He I 5875 A indicates that its line pro- 
file is affected by processes in addition to the surface motions 
produced by the tidal forces. 

A clue to the source of additional variability can be found 
by comparing the ratio of the perturbed and unperturbed theo- 
retical profiles with the observational analogue. The left-hand 
panel of Figure 10 presents the theoretical residuals in the 
phase interval 0.96-1.43, showing the emission "bumps" and 
excess absorption features as they move accross the profile. 
These features are present only within the velocity range ±100 
km s"'. The ratios of observed spectra were constructed using 
the 2001 Jan REOSC ((f) =0.304) spectrum in the denominator 
as a template, due to its apparently normal shape (see Figure 1). 
Before constructing the ratio, the template was shifted in veloc- 
ity by amounts corresponding to the RV of the spectrum to be 
used in the numerator. Examples are shown in Figure 10 (right) 
and compared with theoretical residual spectra at the corre- 
sponding orbital phase. The observational residuals clearly ex- 
tend over a broader velocity range and indicate the presence 
of actual superposed emission, not just the effects caused by 



' The qualitative nature of the variations for a circular orbit is very 
similar. 



ratio-ing line profiles affected by tidal oscillations. This is not 
surprising since we know that Ha is in emission and that its 
line profile is variable (Margon & Bradt 1977; Koenigsberger 
et al. 2003). It is unfortunate, however, that our high resolu- 
tion data do not include Ha preventing us from searching for a 
connection between the Ho- emission and this deduced super- 
posed emission. But a common source for both emissions is 
likely since the Gaussian full-width at half maximum intensity 
of the Ha emission in the 2001 Oct spectrum is 280 km s~', 
which is similar to the width of some of the residual He I 5857 
A profiles. 

Whether the emission is produced in the vicinity of the col- 
lapsed companion or is associated with the B-supergiant can 
be answered with the radial velocity variations, as displayed 
in Figure 1 1 . Although there is a trend for more negative RVs 
around (p =0, just a few of the observations follow this trend, 
and the amplitude is only ~70 km s"'. This argues against the 
residual emission arising primarily in the vicinity of the col- 
lapsed object where large-amplitude RV variations associated 
with its orbital motion would be expected. Hence, we conclude 
that the wind close to the B-supergiant is the most likely source 
for the bulk of the superposed emission. The atmospheres of 
B-supergiants are, in general, very active and the He I 5875 
A line is reported to undergo large RV changes even in single 
stars (Kaufer et al. 2006; Morel et al. 2004), so the superposed 
emission in LSI-i-65 010 is just one more manifestation of the 
activity. 

In summary, the observed line-profile variability in the He 
I 5875 A line exceeds that which predicted from the tidal in- 
teraction model with input parameters as we have adopted for 
LSI+65 010. This is caused by superposed emission which we 
find to be associated primarily with the B-supergiant. The pres- 
ence and variability of this emission indicates significant sur- 
face activity which, however, does not appear to correlate with 
orbital phase and may thus be similar to intrinsic variability that 
is observed for other B-type supergiants (Morel et al. 2004). 

4.3. The 30.7 day "superorbital" period 

If the "superorbital" period discovered by Farrell et al. is re- 
lated to oscillations of the primary star, this period should ap- 
pear in the optical data. Beskrovnaya (1988) reports polariza- 
tion variations with a characteristic timescale of -month, sup- 
porting this possibility. Although this variability timescale is 
not apparent in the RV data, we do find indications of its pres- 
ence in the strength of the residual emission. Figure 12 is a plot 
of the residual emission equivalent width measured by integrat- 
ing over the entire emission feature, showing a possible modu- 
lation with P=30.79 days, which is within the +0.1 day uncer- 
tainty given by Farrell et al. The Farrell et al. To corresponds 
to "minimum modulation of the light curve" which, according 
to their Fig. 3, occurs during the transition from minimum to 
maximum X-ray flux. Hence, the X-ray P super "high" state is at 
phase (p super ~0.25. This coincides with the general trend of the 
data in Figure 12. 

The residual emission line profiles are illustrated in Figure 
13 where they are stacked in order of increasing 30.79-day 
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Fig. 8. Montage of the He I 5875 A average (left) line profiles 
stacked in order of increasing orbital phase from bottom to top. 
The right-hand panel contains line profiles of nights for which 
only one or two spectra were acquired. Dotted profiles are the 
theoretical line profiles, rotationally broadened to v sin(/)=96 
km s whose velocity shifts correspond to a Keplerian orbit 
of a 16 Mq star with a 1.7 M© companion. 

phase. There is a trend for broader and more symmetrical emis- 
sion near phase (bottom of the plot), with a narrower central 
peak developing thereafter. At other phases, two peaks are of- 
ten observed, but the relative intensity of the two peaks can 
change from one night to the next. 

5. Conclusions 

Tidal interactions in non-synchronously rotating binary sys- 
tems produce strong perturbations on the external layers of the 
stars, particularly when the orbital separations are relatively 
small and the stellar radii are large. The binary system LSI-i-65 
010=2801 14-1-650 is precisely such a system. The exploratory 
calculations performed for this paper are aimed at establishing 
the plausibility that the X-ray variability patterns of LSI+65 
010=280114+650 may be attributable to these tidal interac- 
tions. We show that the oscillations on the surface of the B- 
supergiant that are predicted in our calculations have timescales 
that are consistent with the periodic timescales found in X- 
ray data. We propose that these oscillations may modulate the 
B-star's wind, leading to a fluctuating accretion rate onto the 
neutron star, and hence. X-ray variability. Figure 14 is a car- 
toon depicting the scenario for oscillations with P~2-3 hours 
which are assumed to produce enhanced-density shells. This 
idea, originally proposed by Finley et al. (1992) was discarded 
when the optical component of 280114+650 was found to be 
a B-supergiant instead of a Cep-type star, since such short 




-200 -100 100 

Radial Velocity (km s"') 



Fig. 9. Theoretical model line-profiles stacked as a function of 
orbital phase for the tidally-perturbed 16 Mq, Ri =25 star in 
the eccentric orbit. The dotted profiles correspond to the unper- 
turbed line profiles. The strongest variability is observed be- 
tween periastron passage {<p =0.00) and <p ~0.30 and may be 
described in terms of a wave pattern that moves from the blue 
to the red wing. 

timescale pulsations are not predicted for single supergiant 
stars. We rescue Finley et al.'s hypothesis by proposing that the 
presence of a companion can lead to the appearance of such 
short periods. 

The precise mechanisms by which the oscillations may 
translate into a structured stellar wind are not clear, but we 
speculate that an answer may lie in the combination of the 
radial pulsation components and the shear energy dissipation 
due to the tangential motions of the outer stellar layers slid- 
ing over the more slowly rotating inner layers. These two com- 
bined sources (the radial pulsation component and the shear 
energy dissipation), when added to the radiation field that nor- 
mally drives the stellar wind, might cause the inhomogeneities 
to appear on the same timescale as that of the largest amplitude 
oscillations. In order to establish a connection between the tidal 
oscillations and the modulation of mass-ejection in L8I+65 
010, it is necessary to analyze spectral lines that are less suscep- 
tible to contamination from emission, together with observa- 
tions of Ha. It is important to note that such observations have 
been successfully performed in other stars (Fullerton, Gies & 
Bolton 1996; Prinja et al. 2004; Kaufer et al. 2006) showing 
a link between photospheric activity and perturbations in the 
dense inner regions of stellar winds. 

It is interesting to note that the largest oscillation ampli- 
tudes induced by the tidal forces occur near the primary tidal 
bulge. Hence, the periodic modulation of the stellar wind may 
only be present in a small solid angle oriented in the general 
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Table 2. System and parameters for the eccentric (E) and two circular orbit (C1,C2) models. 



Parameter 


Empirical values 


Model-E 


Model-Cl 


Model-C2 


Comments 


P,„-/i fdavs^ 


11.5995±0.0053 


11.5980 


11.5980 


11.5980 


den Hartog et al. 


P„,,„,. (days) 

super \ J J 


30.7 


none 


30.4 


30.7 


Farrell et al. 


Pfi,,^^ ('hours) 

' flare v*^'^"'-''/ 


-2.7 


1.6-3.4 


-2.3 


-3.2 


Hall et al. 


e 


0.0 or 0.16 


0.16 


0.00 


0.00 


Crampton et al. 


T,,, (MJD) 


44134.4 








Crampton et al. eccentric orbit 


Trirr (MJD) 


44133.8 








Crampton et al. circular orbit 


' A— orb v^*-"-""-^/ 


51825.2 








den Hartog et al. "X-ray maximum" 


Tn (MJD) 


43669.9 








Tr(rr"40 cyclcs; elongation 


Ty ,„„„ (MJD) 

A— super v-^'-"" / 


50108.2 








Farrell et al. "minimum modulation" 


Ml (Mo) 


16.0±5 


16 


16 


16 


Reig et al. 


M, (Mo) 


1.7 


1.7 


1.7 


1.7 


Reig et al. 


Ri (Rg) 


37±15 


25 


40 


36.2 


Reig et al. 


V sin i (km s"') 


96±20 


100 


96 


96 


Reig et al. 


i (deg) 


61-90 


90 


65 


90 


Reig et al. 


Yrot (km s"') 


110 


100 


106 


96 




flp^,- (deg) 




27 






This paper 




0.4-0.9 


0.66 


0.615 


0.615 




v(R2 day-') 




0.087 


0.350 


0.120 


viscosity 


5R/Ri 




0.0067 


0.0050 


0.0050 


thickness of external layer 




Fig. 10. Left panel: Ratios of theoretical line profiles of the 
previous figure, illustrating the predicted differences between 
the perturbed and corresponding unperturbed photospheric ab- 
sorption. Orbital phases are indicated. Right panel: Ratios of 
the observed profiles compared to the same theoretical varia- 
tions (dashes) that are illustrated in the right-hand panel. The 
much broader velocity range in the observations (indicated by 
the vertical lines) implies the presence of residual emission, in 
addition to variability of the underlying photospheric absorp- 
tion. 



direction of the neutron star. The localized nature of the oscil- 
lations would explain the lack of a clear ~2.7 hr periodicity 



Fig. 11. Centroid of the residual emission as a function of or- 
bital phase. Stars represent the 1993 data; filled -in squares, the 
1995 data; and open squares all other data except the 2004 Oct 
values which are indicated with crosses within open squares. 
A trend for more negative velocities around (/> -Q is present, 
although the small amplitude suggests an origin near the B-star 
rather than near the collapsed companion. 



in the He I 5875 A data, since the irregular variations that are 
present over the whole B-star surface most likely drown out the 
small-scale periodic effects. Thus, an interesting conclusion is 
that the X-ray emission is a much finer diagnostic of the local 
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Fig. 12. Equivalent width of the residual emission plotted as a Fig. 13. He I 5875 A residual emission stacked in order of in- 



function of the P super -30.79 day phase. A possible modulation 
is present, illustrated with a sine curve. Symbols are the same 
as in Figure 1 1 . Here phase 0.0 corresponds to the time of min- 
imum modulation in the X-ray light curve of Farrell et al., with 



X-ray "high" occuring near i 



=0.25. 



creasing Psuper =30. 79-day phase (from bottom to top), and 
with separations proportional to the 4' super- As in Figure 12, 
zero phase corresponds to "minimum modulation" or the X- 
ray light curve (Farrell et al.). The arrows indicate our deduced 
approximate times of X-ray minimum {(psuper ~0.6-0.8) and X- 
ray maximum {<psuper ~0. 1-0.3). Dots are used to differentiate 
between overlapping spectra on the plot. 



wind structure than the optical line profiles, the latter reflecting 
the integral of the absorption/emission over the stellar disk. 

The tidal interaction model predicts a large difference be- 
tween the pulsation energy dissipation rates, E, near periastron 
and apastron for an eccentric orbit. This implies a significant 
variation in surface activity as a function of orbital phase, as- 
suming that the surface activity is associated with E. X-ray 
variability on the orbital timescale is indeed present, but the 
line profiles in our optical data do not display strong changes 
between periastron and apastron phases, weakening the case 
for an eccentric orbit. In a circular orbit, E remains constant 
over orbital phase and thus no modulation of the stellar wind 
with Porb is expected. If the orbit of LSl+65 010 is circular, the 
observed X-ray variability is Ukely due to eclipse effects, rather 
than to differences in the accretion rates onto the companion. 
A circular orbit would also be consistent with the existence of 
a "superorbital" period. 

Our final consideration concerns the possibility that part of 
the X-ray emission may actually originate in the B-supergiant. 
In particular, the tangential component of the perturbation pro- 
duces a differentially rotating structure in the surface layers that 
one might speculate could lead to magnetic field generation. 
Admittedly, the timescales over which a particular region of the 
star develops this structure may be too short for magnetic field 
amplification mechanisms for massive stars, such as proposed 
by Spruit (2002), to operate. However, it would be interesting 
to explore the mechanisms that could convert tidal energy into 
magnetic field energy, as has been done for low-mass stars (for 
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Fig. 14. Schematic illustration of the scenario for the -2.7 hour 
X-ray periodicity, under the assumption that the oscillations 
near the tidal bulge are responsible for modulating the wind 
of the B-supergiant. 
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example, Zaqarashvili et al. 2002). Furthermore, observational 
evidence for surface magnetic field structures in OB stars now 
exists (Henrichs, Neiner, Geers 2003), so it would be interest- 
ing also to explore the effects of the tidal oscillations on as- 
sumed pre-existing magnetic field structures and determine the 
X-ray emitting characteristics due to these perturbations. 
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